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Background: Previously we demonstrated that DNA can 
act as an enzyme in the Pb 2 *-dependent cleavage of an 
RNA phosphoester. This is a facile reaction, with an 
uncatalyzed rate for a typical RNA phosphoester of 
MO" 4 min" 1 in the presence of 1 mM Pb(OAc) 2 at pH 7.0 
and 23°C.The Mg 2+ -dependent reaction is more difficult, 
with an uncatalyzed rate of MO" 7 min” 1 under compara- 
ble conditions. Mg 2+ -dependent cleavage has special rele- 
vance to biology because it is compatible with intracellular 
conditions. Using in vitro selection, we sought to develop a 
family of phosphoester-cleaving DNA enzymes that 
operate in the presence of various divalent metals, focusing 
particularly on the Mg 2 * -dependent reaction. 

Results: We generated a population of >10 13 DNAs con- 
taining 40 random nucleotides and carried out repeated 
rounds of selective amplification, enriching for molecules 


that cleave a target RNA phosphoester in the presence of 
1 mM Mg 2 *, Mn 2 *, Zn 2+ or Pb 2+ . Examination of indi- 
vidual clones from the Mg 2 * lineage after the sixth round 
revealed a catalytic motif comprised of a three-stem junc- 
tion. This motif was partially randomized and subjected to 
seven additional rounds of selective amplification, yielding 
catalysts with a rate of 0.01 min" 1 . The optimized DNA 
catalyst was divided into separate substrate and enzyme 
domains and shown to have a similar level of activity 
under multiple turnover conditions. 

Conclusions: We have generated a Mg 2 *-dependent 
DNA enzyme that cleaves a target RNA phosphoester 
with a catalytic rate M0 5 -told greater than that of the 
uncatalyzed reaction. This activity is compatible with 
intracellular conditions, raising the possibility that DNA 
enzymes might be made to operate in vivo. 
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Introduction 

The isolation, by in vitro selection [1—3], of a family of 
DNA enzymes that exhibit Pb 2+ -dependent RNA phos- 
phoesterase activity [4], confirmed earlier speculation [5] 
that DNA could function like protein- and RNA-based 
enzymes. More recently, in vitro selection was used to 
isolate a family of DNA enzymes that catalyze the joining 
of two DNA oligomers, one bearing a 5’ hydroxyl and 
the other a 3 f imidazole-phosphate [6], suggesting that 
DNA may be capable of an extensive range of catalytic 
activities. Complex secondary and tertiary structures can 
be formed by single-stranded DNA [7—11] and its 
analogs [12-14]. In addition, chimeric RNA-DNA mol- 
ecules can form active hammerhead ribozymes [15,16], 
and an example exists of a Cu 2+ -dependent self-cleaving 
DNA [17]. Here we report the isolation of new RNA 
phosphoester-cleaving DNA enzymes that function with 
an expanded range of divalent metal ion cofactors, and 
we describe the characteristics of a catalytic DNA that 
operates under physiological conditions. 

A number of different self-cleaving ribozyme motifs 
[18] have been isolated from natural sources and have 
been studied for their potential use as sequence-specific 
‘catalytic antisense’ agents [19,20]. Each requires Mg 2 * 
or some other divalent metal cation as a cofactor to 
accelerate the cleavage of a particular RNA phospho- 
ester. We applied in vitro selection to large pools of 


random-sequence DNA, screening for molecules that 
catalyze this same reaction in the presence of various 
divalent metals. We focused mainly on the Mg 2 *- 
dependent reaction because of its biochemical interest 
and potential relevance to the cellular environment. The 
pK a of a Mg 2+ -bound water is 1 1 .4, making this a chal- 
lenging reaction at neutral pH. A Mg 2+ -dependent 
DNA enzyme that operates in vivo would be required 
to function in the presence of 0.5-3.5mM free Mg 2 * 
[21]. Accordingly, DNAs were selected for their ability 
to cleave a target RNA phosphoester at pH 7.0 in the 
presence of 1 mM Mg 2 *. 

Results and discussion 

A population of >10 13 different double-stranded DNAs, 
each containing a 5-biotin and a single embedded 
ribonucleotide (Fig. la), was immobilized on a streptavidin 
column and subsequently denatured to produce a matrix 
that displays single-stranded molecules. These molecules 
also contained a 40-nucleotide random-sequence domain, 
flanked by two pairing regions that were intended to 
function as substrate-binding domains. RNA-cleaving 
DNAs were isolated by eluting the column with buffered 
solutions containing 1 mM of either Mg 2 *, Mn 2+ , 
Zn 2 * or Pb 2 *, in sequential order. The selected DNAs 
were amplified separately by the polymerase chain reac- 
tion (PCR) and the 5 '-biotin and embedded ribonucleo- 
tide were reintroduced by additional PCR amplification 
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Fig. 1. Isolation of DNAs that catalyze the metal-dependent cleavage of an RNA phosphoester. (a) Starting pool of random-sequence 
DNAs, engineered to contain two substrate-binding domains. Each member of the pool contains a S'-terminal biotin (encircled B), a 
single embedded ribonucleotide (rA) and a 40-nucleotide random-sequence domain (N 40 ). (b) Selective-amplification scheme for iso- 
lation of DNAs that catalyze the cofactor-dependent cleavage of an RNA phosphoester. (c) Self-cleavage activity of the starting pool of 
DNA (0) and populations obtained after either the fifth (5) or sixth (6) rounds of selective amplification, for each of the four metal- 
dependent lineages. [ 5 ' - 3 2 P] - 1 a be I ed precursor (Pre) was incubated at 23 'C for 1 h in buffer B with either none (-) or 1 mM (+) added 
divalent metal ion, as indicated. M is [5'- 32 P] -labeled primer 3, corresponding in sequence to the expected 5 1 cleavage product (Civ). 


(Fig. lb). The entire selective-amplification procedure was 
carried out repeatedly for each of the four lineages. 

After five (Mn 2+ and Pb 2+ ) or six (Mg 2 * and Zn 2+ ) 
rounds, all four lineages showed catalytic activity in the 
presence of their respective metal cofactor (Fig. lc). 
Product formation under the selection conditions for 
the Mg 2 *, Mn 2+ , Zn 2+ and Pb 2 * populations was 3 %, 
11 %, 24 % and 62 %, respectively (Table 1). These values 
are inversely proportional to the pK a values of the cor- 
responding metal-bound water (11.4, 10.6, 9.0 and 7.7, 
respectively). A similar relationship has been seen when 
substituting other divalent metals for Mg 2+ with the 
hammerhead ribozyme [22], in both cases suggesting 
that the metal hydroxide is directly involved in catalysis. 

Examination of 25 clones [23] from each lineage re- 
vealed that the engineered substrate pairing incorporated 


into the original pool was frequently disrupted by muta- 
tions that accumulated during the selective-amplification 
process. Many individuals from each of the selected pop- 
ulations have sequence similarity to the Pb 2+ -dependent 


Table 1. Activities of four catalytic DNA lineages. 


Lineage 


Metal 3 



Mg 2+ 

Mn 2+ 

Zn 2+ 

Pb 2+ 

Mg 2+ 

3 

32 

50 

51 

Mn 2 * 

0 

11 

21 

49 

Zn 2+ 

0 

11 

24 

42 

Pb 2+ 

0 

2 

5 

62 


a Shows % of population cleaved when lOnM precursor was 
incubated for 1 h with 1 mM metal ion. Conditions as for Fig. 1 c. 
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Fig. 2. Mg 2+ -dependent DNA catalysts, (a) Sequence and pre- 
dicted secondary structure of the most active catalyst isolated from 
the Mg 2+ -dependent lineage after six rounds of selective amplifica- 
tion. Underlined sequences correspond to the engineered sub- 
strate-binding domains included in the original pool design (Fig. 
la). Conserved residues are shown in green. Boxed hairpin depicts 
the optimized structure of a representative clone obtained after re- 
selection. (b) Sequence alignment of individuals isolated after ran- 
domization of the most active Mg 2+ -dependent clone and seven 
additional rounds of selective amplification. The sequence of each 
clone matches the prototype sequence except where noted. The 
arrowhead identifies the clone shown in the inset in (a). 

DNA enzyme described previously [4]. The clones iso- 
lated from the Mg 2+ lineage, however, have a more 
complex structure composed of a three-stem junction. 
All of these clones have very similar sequences, with only 
four being unique. The most active clone, which 
occurred five times, is shown in Figure 2a. Under the 
selection conditions, this individual has a first-order reac- 
tion rate, k obs , of 0.002 min'" 1 . It is 3-, 19- and 5-fold 
more active in the presence of 1 mM Mn 2+ , Zn 2+ and 
Pb 2+ , respectively. 

We prepared a degenerate pool of synthetic DNA based 
on the sequence of this most-active clone, randomizing 
the 40 nucleotides that lie between the engineered 
pairing domains such that the original nucleotide 
occurred with a probability of 0.85 and each of the 
other three nucleotides occurred with a probability of 
0.05. We carried out seven additional rounds of selec- 
tive amplification, based on activity in the presence of 
1 mM Mg 2 *, and examined 30 clones from the result- 
ing population (Fig. 2b). This revealed one highly vari- 
able and two strictly conserved regions within the 




Fig. 3. Characterization of the intramolecular reaction catalyzed 
by the optimized catalyst. All reaction mixtures contained 
10 nM [5 ’- 32 P] -labeled precursor and were incubated at 23 ’C 
in 50 mM HEPES (pH 7.0), unless otherwise indicated, (a) 
Effect of NaCI concentration on k bs in the presence of 1 mM 
MgCU. The inset displays a detail of the data for 0-0.25 M NaCI. 

(b) Effect of MgCI 2 concentration on k obs in the the presence of 
1 M (•) or no (O) NaCI. (c) Temperature dependence of k obs in 
the presence of 10 mM MgCl 2 and 1 M NaCI. Precursor DNA 
was prepared by either 5 cycles (for (a) or (b)) or 20 cycles (for 

(c) ) of PCR amplification (see Materials and methods). 
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Fig. 4. (a) Sequence and secondary structure of complexes formed 
between various substrates and enzymes. Conserved residues are 
shown in green, (b) Eadie-Hofstee plot for the reaction involving S2 
and E6. Initial rates of cleavage were determined for reaction mix- 
tures containing 100 nM E6 and either 4, 8, 16, 32 or 64 jxM S2. 

40-nucleotide randomized domain. The sequence varia- 
tion at positions 10-29 is consistent with the existence 
of a stem-loop structure. 


One of the optimized DNA catalysts isolated from the 
reselected pool (Fig. 2a, inset) has a k obs of 0.01 min' 1 
in the presence of 1 mM Mg 24 and 0.02 min -1 in the 
presence of 1 mM Pb 2 *. Addition of 1-10 piM Pb 2 * to a 
reaction mixture containing 1 mM Mg 2 * resulted in no 
significant increase in the catalyzed rate of phosphoester 
cleavage. The rate of the Mg 2 *-dependent reaction is 
10 5 -fold enhanced compared to the uncatalyzed rate 
of cleavage of 5-GTAGAGAAGG rA TATCACT, 
measured under the same reaction conditions. This 
rate enhancement approaches that observed for the 
hammerhead ribozyme [24]. 

At low concentrations, NaCl inhibits the rate of cataly- 
sis, but the rate increases linearly with increasing con- 
centrations of NaCl above 0.1 M (Fig. 3a). This 
inhibitory effect is not detected when lOmM Mg 2 * is 
used. A plot of activity versus Mg 2 * concentration 
reveals a K M for Mg 2 * of lOmM, with a maximum cat- 
alytic rate of 0.08 min" 1 reached at saturation (Fig. 3b). 
The rate of RNA phosphoester cleavage is largely inde- 
pendent of temperature over the range of 15—35 °C, but 
drops sharply at higher temperatures (Fig. 3c). The DNA 
catalyst has a broad pH optimum, with k obs ~10 -2 min' 1 
between pH 6.2 and 8.6 in the presence of lOmM 
MgCl? and 1 M NaCl. Under simulated physiological 
conditions (2mM MgCl 2 , 150mM KC1, 50 mM HEPES 
(pH 7.4), 37 °C) the rate of RNA phosphoester cleavage 
is 0.001 min -1 . 

The optimized DNA catalyst was divided into separate 
substrate and enzyme domains, which were tested for 
activity in an intermolecular context (Fig. 4a). The 
oligomers SI and El correspond precisely to a truncated 
version of the optimized DNA catalyst; E2 contains a 
single A to T change that repairs a putative A*A mis- 
match within one of the two substrate-binding domains. 
Both El and E2 catalyze the cleavage of SI , yielding a 5‘ 
product with a phosphate at its 3' terminus and a 3’ 
product with a free 5' hydroxyl. Under multiple turnover 
conditions, employing 20 fxM SI and 0.2 p,M E2, 
3.4 piM products were formed after a 5 h incubation in 
the presence of lOmM MgCl 2 and 1 M NaCl at pH 7.0 
and 23 °C, corresponding to 17 turnovers of the enzyme. 
This rate closely approximates k obs for the unimolecular 
self-cleavage reaction. 

Constructs E3 and E4 (Fig. 4a) showed no catalytic 
activity, suggesting that mutations within either of the 
two strictly-conserved unpaired regions are not toler- 
ated. Truncation of the substrate (S2),and either short- 
ening or altering the sequence of the central hairpin of 
the enzyme (E5, E6), does not significantly alter enzy- 
matic activity. Cleavage of S2 by E6, in the presence of 
10 mM MgCl 9 and 1 M NaCl, proceeds with a k cat of 
0.039 min -1 and K M of 13 |±M (Fig. 4b). E6 is unable 
to cleave either an all-DNA or all-RNA analog of S2. 

The Mg 2 *-dependent DNA enzyme developed here is 
structurally more complex than the Pb 2 *-dependent 
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molecule generated previously [4]. The latter compound 
is comprised of two stem structures that bind the 
substrate, surrounding a central core of 15 nucleotides. 
Six of the core nucleotides are highly conserved. The 
Mg 2+ -dependent compound also contains two stems 
that are involved in substrate binding and has a central 
core containing 15 highly-conserved nucleotides and a 
stem-loop structure of variable sequence. The secondary 
structure of the Mg 2+ -dependent DNA enzyme is 
similar to that of the hammerhead ribozyme, which also 
catalyzes the Mg 2+ -dependent cleavage of an RNA 
phosphoester. There is no reason to believe, however, 
that the two molecules adopt similar tertiary structures. 

The Mg 2+ -dependent DNA enzyme provides a rate 
enhancement of~10 5 compared to the uncatalyzed reac- 
tion. This is similar to the rate enhancement observed for 
the Pb 2+ -dependent DNA enzyme [4] and is about 
10-fold less than that of the hammerhead ribozyme 
operating in the presence of Mg 2+ [25]. The catalytic 
activity of both the Mg 2+ -dependent DNA enzyme and 
hammerhead ribozyme is significantly lower at physio- 
logical MgCl? concentrations (0.5-3.5mM) than at a 
MgCl ? concentration of 10 mM (Fig. 3b) [22]. We have 
made no attempt to improve the Mg 2+ -binding affinity 
of the DNA enzyme, for example, by carrying out in 
vitro selection at very low MgCl 7 concentrations. Such 
optimization is best deferred until DNA enzymes are 
developed that cleave a target site within a segment of 
RNA that is of biological interest. 


Significance 

We have shown that DNA enzymes can use a 
variety of metal ion cofactors and can catalyze a 
biologically-relevant reaction under conditions 
similar to those that exist within cells. It is now 
clear that the occurrence of DNA enzymes in 
living systems, either extant or extinct, cannot be 
ruled out based on any perceived physical and 
chemical inadequacies of DNA itself. Moreover, 
it may now be practical to tailor-make, through 
in vitro selection, RNA-cleaving DNA enzymes 
that specifically target cellular RNAs and 
that can be used in a therapeutic capacity as an 
alternative to ribozymes. 


Materials and methods 

Oligonucleotides and oligonucleotide analogs 
Synthetic DNAs and DNA analogs were purchased from 
Operon Technologies. Primer 3a, [5'- 32 P]-GGGACGAATT- 
CTAATACGACTCACTATrA-3’, was prepared by 5*-labeling 
primer 3 using [y- 3 -P]ATP andT4 polynucleotide kinase. SI 
was prepared as described previously [4] and S2 was purchased 
from Oligos Etc. 

DNA pool preparation and in vitro selection 

The starting pool of DNA was generated by extending 150 

pmoles of primer 3b, 5'-biotin-GGGACGAATTCTAATA- 


CGACTCACTATrA-3', in the presence of an equal amount 
of template, 5 '-GTGCCAAGCTTACCGTCAC-N 40 - 
GAGATGTCGCCATCTCTTCCTATAGGAGTCGT- 
ATTAG-3’ (N indicates equal representation of G, A,T and 
C), in a 300-jxl reaction mixture containing 0.2 U jxl" 1 Taq 
DNA polymerase, 1.5 mM MgCL>, 50 mM KC1, lOmMTris- 
HC1 (pH 8.3 at 23 °C), 0.01% gelatin and 0.2 mM of each 
dNTP, for 1 min at 92 °C, 1 min at 50 °C and lOmin at 
72 °C.The extended DNA was recovered [4], resuspended in 
buffer A (1 M NaCl, 50mM HEPES (pH7.0), 0.02mM 
EDTA) and immobilized by passing the solution through 
three successive streptavidin affinity columns. The columns 
were washed with five 100-pl volumes of buffer A and five 
100-pd volumes of 0.2 N NaOH, equilibrated with buffer A, 
eluted over the course of 1 h with three 20-pd aliquots of 
buffer B (1M NaCl, 50 mM HEPES (pH 7.0)) with added 1 
mM MgCl 2 and again washed with five 100-p.l volumes of 
buffer A. The process was repeated, successively replacing the 
MgCU with 1 mM MnCl 2 , ZnCl 2 and Pb(OAc) 2 , thereby 
creating four separate DNA lineages. DNA was recovered 
from each eluate by precipitation with ethanol following 
addition of 20 pmoles each of primer 1, S'-GTGCCAAGCT- 
TACCG-3', and primer 2, 5-CTGCAGAATTCTAATAC- 
GACTCACTATAGGAAGAGATGGCGAC-3'. Each of the 
four recovered populations of DNA was amplified separately 
by PCR in a 100-|xl volume over 25 temperature cycles [4]. 
lOp.1 of the resulting mixture was used to initiate a 100-p.l 
nested PCR containing 50 pmoles each of primer 1 and 
primer 3b, which was amplified over six temperature cycles. 
Selective amplification was carried out for five successive 
rounds with Mn 2+ and Pb 2+ , with reaction times of 1 h 
(rounds 1-3), 20min (round 4) and 1 min (round 5), and for 
six rounds with Mg 2+ and Zn^ + , with reaction times of 1 h 
(rounds 1-4), 20 min (round 5) and 1 min (round 6). Rounds 
1—3 of the reselection involved reaction in buffer B with 
added 1 mM MgCl 2 for 1 h. Rounds 4-7 involved pre- 
reacting with buffer B with added 1 mM Pb(OAc) 2 for 1 hr, 
followed by reaction with buffer B with added 1 mM MgCl 2 
for 1 h (round 4), 20 min (rounds 5 and 6) or 1 min (round 7). 

Preparation of precursor DNAs 

Precursor molecules, from either a population of DNAs or 
an individual clone, were prepared with a 5’- 32 P label by 
PCR amplification in a 25-jxl reaction mixture containing 
10 pmoles each of primer lb (5'-biotinylated form of 
primer 1) and primer 3a, 0.5 pmoles input DNA and 0.1 U 
jjlD 1 Taq polymerase, which was incubated for either 5 or 20 
cycles of 1 min at 92 °C, 1 min at 50 °C and 1 min at 72 °C 
under the conditions described above. Precursor DNA pre- 
pared by 20 cycles of PCR amplification consistently showed 
4—5-fold lower activity than that prepared by 5 cycles of 
amplification. Following PCR, NaCl was added to a final 
concentration of 1 M and the entire mixture was passed 
through a streptavidin column. The column was washed 
with buffer A and the single-stranded precursor DNA was 
recovered by elution with 40 pd of 0.1 N NaOH, neutralized 
with an excess of sodium acetate (pH 5.2), precipitated 
with ethanol and purified by denaturing polyacrylamide 
gel electrophoresis. 
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